MOMENTS OF PRODUCTS OF ELLIPTIC INTEGRALS 
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Abstract. We consider the moments of products of complete elliptic integrals 
of the first and second kinds. In particular, we derive new results using elemen- 
tary means, aided by computer experimentation and a theorem of W. Zudilin. 
Diverse related evaluations, and two conjectures, are also given. 
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1. Motivation and general approach 

We study the complete elliptic integral of the first kind, K{x), and the second 
kind, E{x), defined by: 



Definition 1. 
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As usual, K\x), E'{x) are used to denote the function value at the complemen- 
tary modulus, x' = \f\ — x^. Recall that ^F^ denotes the generalized hypergeomet- 
ric series, 
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The complete elliptic integrals, apart from their theoretical importance in arbi- 
trary precision numerical computations ([8]) and the theory of theta functions, are 
also of significant interest in applied fields such as electrodynamics ([16]), statisti- 
cal mechanics, and random walks. Indeed, they were first used to provide explicit 
solutions to the perimeter of an ellipse (among other curves) as well as the (exact) 
period of an ideal pendulum. 

The author was first drawn to the study of integral of products of K and E in 
[9], in which it is shown that 

2! K{xfdx= [ K'{xfdx, (3) 
Jo Jo 

by relating both sides to a moment of the distance from the origin in a four step 
uniform random walk on the plane. 

A much less recondite proof was only found later: set x = (1 — t)/(l + 1) on the 
left hand side of (3), and apply the quadratic transform (4) below, and the result 
readily follows. 
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The four quadratic transforms ([8]), which we will use over and over again, are: 

= T^k(\^) (4) 



1 + X \ 1 + X 

E'{x) = (1 + x)E - xK'{x) (6) 

In the following sections we will consider definite integrals involving products of 
K, E, K', E', especially the moments of the products, with the intent of producing 
closed forms whenever possible. 

The somewhat rich and unexpected results are perhaps surprisingly easy to 
discover, thanks to methods in experimental mathematics, for instance, the integer 
relations algorithm PSLQ, the Inverse Symbolic Calculator (ISC, now hosted at 
CARMA, [12]), the Online Encyclopedia of Integer Sequences (OEIS, [15]), the 
Maple package gfun, Gosper's algorithm (which finds closed forms for indefinite 
sums of hypergeometric terms, [13]), and Sister Celine's method ([13]). Indeed, 
large scale computer experiments ([2]) revealed that there was a huge number of 
identities in the flavour of (3). Once discovered, many results can be routinely 
established by the following elementary techniques: 

(1) Connections with and transforms of hypergeometric and Meijer G-functions 
([16]), as in the case of random walk integrals (Section 4). 

(2) Interchange order of summation and integration (Section 3). 

(3) Evaluate a closed form for the primitive (Section 5). 

(4) Change the variable x to x' (Section 5). 

(5) Change the variable followed by quadratic transforms (Section 5). 

(6) Use a Fourier series (Section 6). 

(7) Apply Legendre's relation (Section 7). 

(8) Differentiate a product of elliptic integrals and integrate by parts (Section 
8). 

2. One elliptic integral 

The moments of K, E, K' , E' are established, for instance, in [S]. Moreover, in 
[6], all moments of the generalized elliptic integrals are similarly derived. Here we 
slightly generalize the results in [8]. 

By a straight-forward application of the formulae (see [6]) 
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These also encapsulate the moments for K\ E', because hj x ^ x', we have 



x"'x"'K'(x) dx 



X 



dx, / xVE'(x)dx 



We note that for convergence, m > — l,n > —2. When m = 1 both formulae 
reduce to a 2-F1 and can be summed by Gauss' theorem ([!]). 

If in addition 2m + n + 1 = in (8), then Dixon's theorem ([!]) applies and 
we may sum the 3F2 explicitly in terms of the F function. For instance, we may 
compute Jq K{x)/x' dx (which also follows from the Fourier series in Section 6). 

In (9), Dixon's theorem may only be applied to the single special case x'E(x) dx. 

3. One elliptic integral and one complementary elliptic integral 

Here we take advantage of the closed form for moments of K', E' which follow 
from (8) and (9): 
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We also use the (hypergeometric) series for K, E which follow from the defini- 
tions: 

F(A; + 1/2)2 x^'^ 
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Hence, for the proposition below, we simply write the elliptic integral as a sum, 
and observe that the order of summation and integration may be interchanged. 

Proposition 1. We have the following moments: 
(1) 
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In this case, odd n yields a closed form in terms of vr'^, and possibly a 
rational multiple of vr, as we can expand the summand (which is a rational 
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function) as a partial fraction much like in Remark 2. Amazingly, all the 
rational multiples of it vanish as we show in Lemma 1. 

(2) 

'1 7r2r(i(n + l))2 /_i 1 n+l n±l 
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In this case, odd n yields a closed form in terms of vr and vr'^. From 
Lemma 1 and (54), the constant of the vr term is ^^^^-^ ■ 

(3) 

In this case, odd n yields a closed form in terms of vr and vr'^. From 
Lemma 1 and (54), the constant of the vr term is 

(4) 

X ^ ^ 8 (n + 2)r(l(n + 2))2^'^^^ l,n±2,ii±4 ^^^^ 

In this case, odd n yields a closed form in terms of vr and vr'^. By using the 
derivative of (1 — x'^)x'^"'E{x)E'{x) and the linear terms above, we deduce 
that the linear term here is also TrTT\- 

4(n+l) 

We observe that E', K' do not have 'nice' expansions around the origin, and 

K do not have 'nice' moments, so the method above cannot to extended to 
evaluate moments of other products. 

Note that by using the symmetry between parts 2 and 3, as well as by applying 
Legendre's relation (38), we obtain linear identities connecting these aF^''s. 

We now prove the observation in part 1 of Proposition 1. 

Lemma 1. When n is odd, the nth moment of K{x)K\x) is a rational multiple 
o/vr^. 

Proof. We experimentally discovered that, letting g{n) := x'^"~^K{x)K'{x) dx, 
we have 

2n^g{n + 1) - (2n - l)(2n^ - 2n + l)g{n) + 2{n - lfg{n - 1) = 0. 

This (contiguous) relation, once discovered, can be proven by extracting the sum- 
mand, simplifying and summing using Gosper's algorithm. Thus, after computing 
2 starting values, the claim is proven. Note that the recursion also holds when n 
is not an integer. □ 

Experimentally, we found that the sequence h{n) := 7r^lQ'^~^^g{n + 1) matched 
entry A036917 of the Online Encyclopedia of Integer Sequences; indeed, they share 
the same recursion and initial values. Moreover, the OEIS provides that 

^, , ^f2n-2kyf2kY 16"r(n + l/2)2 / -n, -n, i, ^ \ . 



MOMENTS OF PRODUCTS OF ELLIPTIC INTEGRALS 5 

The first equality is routine to check as we can produce a recurrence for the 
binomial sum - for instance, using Sister Celine's method; the second equality is 
notational. 

The generating function for h{n) is simply 

oo . 

Y^hin)r = -K{4Vtf, 

n=0 

which is again easy to prove using the series for K{t). Recall that h{n) is related 
to the moments of K{x)K'{x), and thus we have: 

Theorem 1. 

I' Y^,K{x)K\x) dx = ^-K{tf. (15) 
Equation (15) seems to be a remarkable extension of its (much easier) cousins, 

^'^3^/^'(x)dx = |jat), /r3W^'(^^^^ = l^^^^^')~^('))- 

Note that manipulations of (15, 16) give myriads of integrals, we list some of 
which below (below G denotes Catalan's constant): 

r^^^^K'{x)dx = nG, 
Jo X 

I -K{x)K'{x){K{x)- E{x))dx = [ K{xfE'{x)dx, 
Jo ^ Jo 

j\F^(^'^ x'^xK{x)K\x)dx = ^^^^-^ j\-K{tfdt. 

The last identity specialises to 

C — i^^^^^^i^(x)A"(x) dx = -7rC(3). 
Jo ^ 8 

4. Two COMPLEMENTARY ELLIPTIC INTEGRALS 

In [17], Zudilin's Theorem connects, as a special case, triple integrals of rational 
functions over the unit cube with generalized hypergeometric functions (rFg's). 
We state a restricted form of the theorem which is sufficient for our purposes here: 

Theorem 2 (Zudilin). Given ho, ... ^h^ for which both sides converge, 
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In [9], this theorem is apphed to derive closed form evaluations for moments of 
random walks from their integral representations. 

The idea here is to work backwards, writing a single integral involving products 
of elliptic integrals as a double, then a triple integral of the required form, and 
then apply Theorem 2. To do so, we require the following formulae, which can be 
readily verified (see e.g. [14]): 



a/x(1 - x){a - x) 



a \ \ a 



= 2^aE{^\, (19) 



dy 



, ^= = 2K\^d), (20) 

a ^y{l - y){y - a) 

['^==^==dy = 2E\^a). (21) 

Though the simple cases in this section (corresponding to n = 0) are tabulated 
in [11], the general results appear to be new. 
Using the above relations, we have, for instance. 



1 rl 




^ Axdydz 



(1 — y)z{l — z) V x{l — x) 



y x{l - x){l -y)z{l - z) 



The first equality follows from (21), the second from the change of variable 
^-)■ yz, the third from (19), and the fourth from z ^ 1 — z. Now Theorem 2 

applies to the last integral. This strategy works for the rest of this section, so we 

are led to: 

Proposition 2. For all n where the integral converges, 



\n.p.. .2. _ 2^> + mn + 3)^r(i(n + l))^ f-ilil'^ 
^xi^[X)ax iQ(^ + 2)3(n + 4) T{n + 1)^ i, "±3, n±2^ n± 



n+3 n+3 n+7 
2 ' 4 
n+3 n+2 n+4 n+4 n+6 

'0 xuv/t -r ^; V"- "T ± v" "T -l; 

(22) 

When n is odd, the hypergeometric reduces to known constants only involving 
C(3), which we prove in Theorem 3 below. 

Similarly, by building up the E' integral then K', we obtain: 
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Hence, we have: 

Proposition 3. For all n where the integral converges, 







1 94"Yr) -H^2 p H/'^ I / _l I I I n+l n±l B+b 

^^^ ^ 16(n + 2) r(n+l)4 ^l, Il±l, n±2, B±2^ n±2^ n±4 

(23) 

Alternatively, by building up the K' integral then E', we obtain: 

E'(x)K'(x)dx=- [ [ [ ^ dxdyd^, 

" Jo Jo - x)(l - y)z{l -z){l- x{l - y{l - z))) 

So the next proposition allows us to produce some equalities for the jFq's in- 
volved: 

Proposition 4. For all n where the integral converges, 

"1 2^"(v -^l)^(ri -^'^)T (-(n+ 1))^ /1113ri±3n±3n+7 

x''E'ix)K'{x) dx = ^ I V r .J iF, ( ^: ^: \, 

16(n + 2)^ r(n + l)4 ^ ^ ^1, ^ 

(24) 

Finally, we also have 

1 1 r r r 1 







K'(xf dx = - I I I , dxdydz, 
" Jo Jo y/xil - x)yil - y)zil - z){l - x{l - 2/(1 - z))) 



Therefore: 

Proposition 5. For all n where the integral converges. 



^nT^I( \2^ ^ ^^2^ F i 2' 2' 2' 2' 2 ' 2 ' 4 

^ ^ 16 T{n + IY ^'^^^l,«±i,B±2,B±2^n±2^n±2 

(25) 

Again, the odd moments are computable in terms of C(3), with the particularly 
simple 

^ xK'ixfdx = ^C(3). (26) 

Remark 1. Therefore, all the odd moments of K'^,E'^,K'E' have particularly 
simple forms. By using (29), we can iteratively obtain all the odd moments of 
K^, E^, KE. As an example, x^Kixf dx = |(2 + 7C(3)). 

Now we prove the observation regarding the appearance of C(3): 
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Theorem 3. When n is odd, the nth moment of K''^ , E''^ , K' E' , K'^ , E'^ and KE 

IS expressible as a + b({3), where a, 6 e Q. 

Proof. We prove the case for the pair K'^ and K^; the other 2 pairs are similar. 
Firstly, when n is odd, the summand of the jFq for K'^ is a rational function: 

{2k + m + l){k + l)\k + 2)^ ■ ■ ■ (fc + mf 
{k + l/2y{k + 3/2)^ ■■■{k + m+ 1/2)^ ' ^ ^ 

for we have ignored any rational constants and wrote n = 2m + l. We can explictly 
sum (27) and verify the statement of the theorem for the first few moment of K'"^. 
By using (29), we can likewise do this for K^. 

We now use the recursion (51) to see that the statement holds for all odd 
moments of K^. By using (29) (but with the role of x and x' reversed), we obtain 
the result for K'^. 

□ 

Remark 2. We sketch another proof by expanding (27) into partial fractions. 

As each partial fraction has at most a 4th power on the denominator, the sum 
can involve at most C(2), C(3), C(4), and some contribution from the linear denom- 
inators. 

Because the linear denominator terms must converge, their sum must eventually 
telescope, and hence contribute only a rational number. 

We recall that partial fractions can be obtained via a derivative process akin to 
computing Taylor series coefficients; indeed, if we write 

fix) _ A, ^ A2 ^ ^ An 

(x — a)" (x — a) (x — ay (x — a)" ' 

then An = f{a), = /'(a)/l!, ...,A, = f^^~'\a)/{n - 1)!. 

When applied to (27), it is easy to check that, when n = 3 mod 4, the presence of 
the numerator 2k+m+l makes the terms with quadratic and quartic denominators 
telescope out, leaving us with rational numbers (these terms occur in pairs related 
by the transformation k 1— )■ (— n — l)/2 — k, where said linear numerator switches 
sign). Similarly, the terms with cubic denominators double. 

When n = 1 mod 4, 2k + m + 1 cancels out with one of the denominators, 
making it a cubic; it can be similarly checked that its partial fraction can have no 
quadratic term: this is equivalent to showing (27) with all powers of 2A; + m + 1 
removed has derivative at = (— n — l)/4, which is true because it is symmetric 
around that point. 

Therefore only the cubic terms remain, giving us C(3)- 

This type of partial fraction argument is at the heart of the result that infinitely 
many odd zeta values are irrational (see [5], which, incidentally, is the motivation 
for Zudilin's Theorem). 

We note in passing that by either one of the two known transformations for 
non-terminating jFq^s ([4]), we can write our iFq as the sum of two 4-F3's with F 
pre-factors, where one term is easily reduced to closed form when n is odd, while 
the harder term becomes obtainable in light of our reduction. 
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5. Sporadic results 

We list some results found by ad-hoc methods; some are not moment evaluations 
while others are preparatory for later sections. 

5.1. Explicit primitives. A small number of integrals have explicit primitives; 
we list some here: 

Ej^) E{x) E{x) Kj^) Kj^) Kj^) 

^^v^)' T^' T±^' T^' ^7^' ~J^' ■ 



We note in passing that even some good CAS may struggle with finding all prim- 
itives for this short list. 

5.2. Complementary modulus. Using the obvious transformation x i— )■ x', we 
have 

1 /■! 

xK{xyE{xfK'{xyE'{xYdx= / xK'{xyE'(xYK{xyE{xYdx, (28) 

Jo 

More generally: 

[ x^''^^K{xYE{xfK'{xyE'{xYdx= I x{l-x''YK'{xYE'{xfK{xYE{xYdx, 
Jo Jo 

(29) 

and for any function /, 

xx' f{x) dx = x^/(x') dx. 



5.3. Imaginary argument. In [i 1] volume III, some integrals with the imaginary 
argument ix are considered, e.g. 



/•I , 1 

/ xK'{x)K{ix) dx = —Gtx. 
Jo 2 



This can be proven by expanding xK{ix) as a series and summing the moments 
of K'{x). Other evaluations are done similarly. For instance, we can easily obtain 
recursions for the moments of K{ix), E{ix). 
We record here that 

E{i^) = Vk+1 E{^/k/{k + l)), and K{iVk) = 1/Vk + l K{^/k/{k + 1)). 
5.4. Quadratic transforms. Using the quadratic transforms (4, 5), we obtain 

K{xYdx = -/ K'{tY{^] dt, 

(30) 



^ K'ixYdx = 2 j K{tY{l + tY~^dt. 



Setting n = 1 we get the known special case 

^ K{x) , vr 



x + l 



2 

dx = — . 
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Using a cubic transform of the Borweins ([7]), this has been generaUzed in [6]. 
The appropriate generahzation of (3) is 

\2 „i /12n2 

dx = 3 / 2F1' '^'^ 



1 /IS 

3' 3 



2 



1 



X dx. 



Equation (3) itself is obtained by setting n = 2 in (30). 
Using (5) on the integrand xK{x)^, we get 

[ 2{1- x)K{xfdx = [ xK{xfdx, 
Jo Jo 

so combined with (30), we deduce 

[ K\xfdx=—! K{xfdx = h! xK{xfdx = h! xK\xfdx. (31) 
Jo 3 Jo Jo Jo 

Using (6, 7), we have 

/•I on+l /■! 

/ E{xY-—-—-^dx= / {E\x) + xK\x)r dx- 
Jo \^ + -Lj Jo 

when n = 1 we can evaulate this in closed form: 

^ Eix) , vr^ 1 

- — "^—L— dx = h -. 

(x + l)3 32 4 

The case n = 2 gives a messier closed form, in view of our knowledge of the 
moments. 

Similarly, it follows from (6, 7) that 

I ^^(^)" (^ + i)n+2 = / (2i?(x)-(l-x2)i^'(x))" dx; 
with the special case 

^ E'(a;) , G 5 

n T = — -I- — 

(x + l)3 8 16' 

5.5. Relation to random walks. In [9], many relations are derived while com- 
puting W4,{n), the nth moment of the distance from the origin of a 4-step uniform 
random walk on the plane. For instance, we have: 

W^il) = ^ / (1 - 3x^)K'{xf dx. 

Jo 

In [10], the following identities are proven, via connections with Meijer G- 
f unctions: 

7]-3 /-VZ r-K/2 

— Wi{-1) = / K{smtYsmtdt = 2 K{smty costdt = / K {sin t)K {cos t) dt, 
4 Jo Jo Jo 

and 



[ K{x)K'{x)x' dx = I K{xfdx. (32) 
Jo Jo 
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6. Fourier series 
As recorded in [3], we have the following Fourier (sine) series valid on (0,27r): 
Lemma 2. 

n=0 ^ ' 

For completeness, we sketch a proof here: 

Proof. By symmetry we see that only the coefficients of sin((2n + l)t) are non-zero. 
Indeed, by a change of variable cos t i— ?■ x, the coefficients are 

4 I' sm((2n+Md, 

The fraction in the integrand is precisely U2n{'x), where Un{x) denotes the Cheby- 
shev polynomial of the second kind ([16]), given by 



k=0 ^ 



k 



We now interchange summation and integration, and use the moments of K'. The 
resulting coefficient contains a 3F2, which after a transformation ([1], section 3.2) 
becomes amenable to Saalschiitz's theorem ([!]), and we obtain (33). □ 

We may apply the same method to obtain a Fourier sine series for E{sint) valid 
on (0,7r), which we had not been able to locate in the literature. In mirroring 
the last step, the resulting 3F2 is reduced to the closed form below using Sister 
Celine's method: 

Lemma 3. 

, , ^r(n + 1/2)2 ^ (n + l/2)r(n+ 1/2)2 

n=0 n=0 

(34) 

Using Parseval's formula on (33) and (34) gives 

7r/2 nTT/2 



/ K{smtydt = 2 K {sin t)E {sin t)dt 



Jo Vi 



X 



K{x)K'{x)dx = -,F,h'^^'^Y 



(35) 



We also get JJ^^^ K{sintY cos(4t) dt as a sum of three 4-F3's, and j^^"^ E{sintY dt 
as a sum of four 4F3's. Section 3.7 of [3] provides more identities with more exotic 
arguments, as well as connections with Meijer G-functions. 
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Experimentally we found that 

-.7r/2 1 rl 



/ -KismtYsmAtdt = / K{xf{x-2x^)dx= xK{x){{l - x'^)K{x) - E{x)) dx 
Jo 4 Jo io 

= [ K'{x)\2x^ -x)dx= [ xK'{x){x'^K'{x) - E'{x))dx 
Jo Jo 

1 

= / xK^{x) - 2xE{x)K{x) dx = — . (36) 
Jo 2 

The last line, an evaluation which is somewhat surprising, is equivalent to 



xK{xy + 2xE{xy - 3xK{x)E{x) dx = 0, 

which is routine to check as we know all the odd moments. 

Inserting a factor of cos(t)^ before squaring the Fourier series (33) and integrat- 
ing, and we are led to 

|.'A-w^d.^^f (23F,(^;y;^|i) -i) ^/'jfjA-W'd.. 

The Fourier series (33) combined with a quadratic transform gives: 

fC sint dt= / , ^ ' dx= / , ^ ' dx = / -^^^dx 
Jo Jo Jo 

' ^'*-''-'d.^A'aV^^rfiV. (37) 



2 Jo IStt 
This result has been generalized in ([6]). 

7. Legendre's relation 

Legendre's relation ([8]) is related to the Wronskian of K and E, and shows 
that the two integrals are closely coupled: 

E{x)K'{x) + E'{x)K{x) - K{x)K'{x) = ^. (38) 

If we take Legendre's relation, multiply both sides by K'{x) and integrate, then 
replacing known moments by their closed forms, we arrive at 

3E'{x)K'{x)K{x) - K{x)K'{xf dx = —. 

8 

Similarly, if we had multiplied by K{x), the result would be 

1 .1 

?>E{x)K{x)K'{x)-2K{xfK'{x) dx = 2E'{x)K{xf-E{x)K{x)K'{x) dx = nG. 

Jo 

(39) 

Using closed form of the moments, we also have: 

/ 2xE'{x)K{xfK\x) - xK{xfK'{xf dx = — , 
Jo 32 
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/ 2xE'{xfK{x)E{x) - xK{x)K'{x)E{x)E\x) dx = — + . 

Jo 16 128 

We can multiply Legendre's relation by any function whose integral vanishes on 

the interval (0, 1) to produce another relation. Suitable candidates for the function 

include x{K{x) - K'{x)), 2{2K'{x) - 3E'{x)), 2E'{x) - K'{x), 2E{x) - K{x) - 1, 

and a vast range of polynomials. For instance one could obtain 

»i 

2E''^{x)K{x) + 2E{x)E\x)K{x) - 5E'{x)K{x)K'{x) + K{x)K'{xYdx = 0. 



Unfortunately, we are not able to uncouple any of the above sums and differences 
to obtain a closed form for the integral of a single product. 



8. Integration by parts 

The following simple but fruitful idea is crucial to this section. We look at the 
rivative {1 — x'^)'^ -^{x'^ K (x)"" E (x^ K' (xY E' (xY) , and integrate by parts to yield 

\l - x^Y-^{x''K{xYE{xYK'{xYE'{xY) dx 
1 

2nx(l - x^Y~^x^K{xYE{xYK'{xYE'{xY dx, (40) 

(where an offset by a constant is possible if the function is a powers of E or E'). 

In practice, we take n, k & {0,1, 2} to produce the cleanest identities. Later, we 
will also explore when n is a half-integer, as well as replacing 1 — by 1 — x in 
(40). 

8.1. Table for products of two elliptic integrals. We now systematically 
analyse the tables kindly provided by D. H. Bailey, the construction of which 
are described in [2]. The tables contain all known (in fact, almost certain all) 
linear relations for integrals of products of up to k elliptic integrals {k < 6) and a 
polynomial in x with degree at most 5. In this subsection we exclusively look at 
when k = 2. 

We use the derivative x-^E{xY = 2E(xY — 2E{x)K{x), and integrate by parts 
to deduce ^ 

/ 'iEi^xf -2E{x)K{x)dx = 1. (41) 
Jo 



More generally, 

1 = [n + k + l] 

Two special cases of the above are prominent: 

»i 

hx^E{xY - 2E{x)K{x) dx = 1, (43) 



[ x^'EixYdx-n [ x''E{xY~^K{x)dx. (42) 
Jo Jo 



and ^ 

/ (n + 2)x""^E(x)2 - 2x'^~^E{x)K{x) dx = 1. (44) 
Jo 
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The derivative of K{x)E{x) (via integration by parts) gives 

;i - ^x^)E{x)K{x) + E{xf - (1 - x^)K{xf dx = 0, (45) 



while more generally 

nx'^-^E{x)K{x)-{n+2)x^+^E{x)K{x)+x^-^E{xf-x^-^K{xf+x''+^K{xf dx = 0. 

(46) 

The derivative of K{xY produces 

/" (1 + x^)K'^{x) dx = 2 [ K{x)E{x) dx, (47) 
Jo Jo 

and more generally, 
»i 

2x''-^E{x)K{x) + {n- 2)x''-^K{xy - nx'^^^K{xf dx = 0. (48) 



The derivative of E'{x)'^ gives 

[ 2xE'{xf - xE' {x)K' {x)dx = [ 2xE{xf - xE{x)K{x) dx = -. 
Jo Jo 2 

The derivative of K'{xY gives 

f 2K'{x)E'{x) - (1 - x^)K\xf dx = 0, 
Jo 

reconfirming a result from random walks ([9]), which was first proven in a much 
more roundabout way via a non-trivial group action on the integrand. 
The derivative of E'{x)K'{x) gives 

[\l - 3x^)E'{x)K'{x) dx = f E'{xf - x^K'{xf dx, 
Jo Jo 

which when combined with the last result, gives 

f {l + ?>x^)E\x)K'{x)dx= [ K'{xf - E'{xfdx. 
Jo Jo 

The derivative of K{x)K'{x) gives 

! x^K{x)K\x) + K{x)E\x) - K\x)E{x) dx = 0, 
Jo 

which, when combined with Legendre's relation (38), results in 



[ 2E'{x)K{x) - (1 - x'^)K{x)K'{x) dx 
Jo 



71 

2' 



Our results above and in the previous sections actually provide direct proofs of 
most entries in the table where the polynomial involved is of degree at most 1. In 
fact, it would simply be a matter of tenacity to prove many other entries in the 
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table where the polynomial is of higher degree. As an example, we indicate how 
to prove one of the few entries that requires more work: 

/ E{x){?>E\x) - K'{x)) dx = -. (49) 
Jo 2 

We note that the left hand side can be written as two 4-F3's. We combine their 
summands into a single term and simplify; the result can be summed explicitly by 
Gosper's algorithm, and the correct limit on the right hand side follows. 

There is only one entry in the tables for which we do not possess a proof: 

Conjecture 1. 

/ 2K{xf - AE{x)K{x) + ?>E{xf - K'{x)E'{x) dx '= 0. (50) 
Jo 

We note that, among moments of products of two elliptic integrals, there are 
only five that we do not have closed forms of: 

E{xf, x^E{xf, E{x)K{x), x^E{x)K{x), x^K{xf, 

as all the odd moments are known, and the other even moments may be obtained 
from these ignition values. Unfortunately, we can only prove four equations con- 
necting them, namely (41, 43, 45, 47). A proof of (50) would give us enough 
information to solve for all five moments. 

8.2. Recurrences for the moments. As already hinted in Lemma 1, the mo- 
ments enjoy recurrences with polynomial coefficients. As an example, by combin- 
ing (44, 46, 48), we may obtain, with Kn ■= x"'K{x)^ dx, 

{n + lfKn+2 - 1n{r? + l)i^„ + (n - l)^i^„_2 = 2. (51) 

This then shows that Kn is a rational number plus a rational multiple of C(3) for 
odd n. 

Similarly, recurrences for other products may be obtained in the same way, 
though the linear algebra becomes more prohibitive. We are able to obtain, as 
another example, for E^ ■= x^E(x)'^ dx, 

{n + l)(n + 3)(n + 5)E„+2 - 2{n^ + + n + 1)E„ + {n - lfEn-2 = 8, (52) 

while the recursion for the moments of EK is straightforward from this using (44). 
The recursion for the moments of K'"^ is identical that (51) except the right hand 
side is 0, not 2. In particular, this approach is used in the proof of Theorem 3. 

8.3. More results. We also discovered some results not found in the tables by 
incorporating constants such as vr and G into the search space. Below we highlight 
some other pretty formulae. Using (1 — x'^)-^{x'^K{xy) gives 



f x^K{xf - xK {x)E{x) dx = 0, 
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while the derivative of x^K{x)K' {x) gives 



I xK{x)K'{x)dx = [ 2x^K{x)K'{x)dx 
Jo Jo 



,„ 1 + 16 
The derivative of x'^K{x)"' imphes 



K{x)K'{x) dx = —. (53) 



1 

x''-^K{x)''-\nE{x) - 2x^K{x)) dx = 0. 
The derivative of x'^^K{x)K'{x) together with Legendre's relation (38) gives 



x^''-\E'(x)K(x) + n(x^ - l)K(x)K'(x)) dx = —. (54) 

8n 



In (40),we can take n = ^. For instance, the derivative of xx'K{x) gives 

f^x^K(x) , 
Jo X 



T 

1 ^2 , 



and the derivative of xx'Kixf' recaptures (35). 
The derivative of x'K[x) gives 

^K{x)-E{x) , vr 

; dx = — , 

XX 2 

note that each part does not converge. 
In fact, 

„2 / 1 3 

x' \ , (55) 



2' 2 
2 



K{x)-E{x) = —2Fr 
therefore 

Jo X 2 Jo a;^ Jo xx' 2 

The derivative of x(l — x)ii'(x)^ gives 

P^^'W^Wd.^ /'ifW^d.. 

io ^ + 1 io 
while that of x(l — x)£'(a;)ii'(a;) gives 

/ §^ " / ^^"^ ~ l)^(x)i^(x) + (x - l)/s:(x)^ dx. 
Collecting what we know about the integral of ii'(x)^, we have the following: 
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Theorem 4. 



1 . rl 

K{xfdx = - / K'ixfdx 

2 In 



C K\xf-dx 
Jo x' 
1 

K{x)K'{x)x' dx 
' 2K{x)E{x) 



x + l 
2 arcsinx 



TT Jq a/T- X 



K{x)K'{x) dx 



1 

arctanh(x)-R'(x)-R''(x) dx. 



Proof. The last two equalities follow from (15); the rest had been proven elsewhere. 

□ 



8.4. Table for products of three elliptic integrals. We now consider the 
linear relations involving the product of three elliptic integrals (/c = 3 in the 
tables). As the number of relations provided by Bailey's tables is huge, we restrict 
most of our attention to a class of integrals that turn out to be mutually related 
by a rational factor. 

Below we tabulate all the products for which neat integrals may be deduced by 
differentiating them and integrating by parts: 



Product: 


Integral: 


K{xf 


jlK{xfE' 


K{xf - m{xfE{x) dx 


= 


K{xfK'{x) 


[x) + K{xfK'{x) - 2K{x)K'{x)E{x) dx 


= 


K'{xfK{x) 


/o E{x)K'{xf - 2E'{x)K{x)K'{x) dx 


= 


K\xf 




K\xf - 3K'{xfE'{x) dx 


= 


E\xf 




/o 5xE'{xf - 3xE'{xYK\x) dx 


= 1 


E{xf 




AE{xf - 3E{xYK{x) dx 


= 1 



We now prove 



[ K{xfK\x)dx = -f K{x)K\xfdx. 
Jo 3 Jq 

We make a change variable x t— )■ {1 — x)/{l + x) to the left hand side, use a 
quadratic transform, then apply x i— )■ 2y/x/(l + x) to part of the result followed 
by another quadratic transform (4, 5), obtaining 

1 .1 

3xK{x)K'{xfdx = / K{x)K\xYdx, 
Jo 

now combining the pieces proves the claim. 
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If we make the change of variable x h-t- (1 — x)/(l + x), then apply (6), we have 



1 + Sio 



Using the derivative of — x)K{x)^, we can show that xK{x)'^E(x)/{x + 
1) dx is also linearly related to the above integral. 

Therefore, gathering the results in this section and equation(31), we have de- 
termined: 

Theorem 5. Any two integrals in each of the following groups are related by a 
rational factor: 

K{xf, K\x)\xK{x)\ xK'{x)\ K{xfE{x), K\xfE\x), ,^'(^)^(^) ^K\x)E{x) 



1 + X 1 + X 

K{x)K'{x)\ K{xfK\x),xK{x)K\xf, xK{xfK\x). (56) 

We cannot yet, however, show that the two groups themselves are related by a 
rational factor, though it is true numerically to extremely high precision. In fact, 
the Inverse Symbolic Calculator gives the remarkable evaluation: 

Conjecture 2. 

Once proven, this also gives explicit closed forms for the integrals of E'{x)K'{x)K{x) 
E{x)K'{x)K{x), and E'{x)K{x)'^, for we can relate each of these to one of the 
above and a constant (by the results of Section 7). 

In view of Theorem 5, equations (10) and (14), interchanging the order of sum- 
mation and integration gives an equivalent form of Conjecture 2: 

i, l,n + 1 



E ° p ( 2' 2 

_j2n + l)2^ l,n + |,n + | 

^ T{n + 1/2^ ( U,-n,-n 



247r4 ^ 



8.5. Products of four elliptic integrals and conclusion. If we take the de- 
rivative for K{x)^, use the integral (30) connecting K'{xY and K{xY, plus a 
quadratic transform, then we obtain 

f 2AE{x)K{xf -%K{xY - K'{xfdx = (59) 

which is one of the first non-trivial identities in Bailey's table for k = 4. Many 
more tabulated relations for products of three and four elliptic integrals can be 
proven, albeit the complexity of the proofs increase. As perceptively noted in [2], 
"[it] seems to be more and more the case as experimental computa- 
tional tools improve, our ability to discover outstrips our ability to 
prove." 
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